arXiv: 1506.09108vl [astro-ph.SR] 30Jun2015 


Resonant Absorption of Transverse Oscillations and Associated 
Heating in a Solar Prominence. II- Numerical aspects 

P. Antolin^, T. J. Okamoto^’®, B. De Pontieu^’^, H. Uitenbroek^, T. Van Doorsselaere®, T. 

Yokoyama^ 

^National Astronomical Observatory of Japan, Osawa, Mitaka, Tokyo 181-8588, Japan 
‘^ISAS/JAXA, Sagamihara, Kanagawa 252-5210, Japan 
^Lockheed Martin Solar and Astrophysics Laboratory, 3251 Hanover Street, Palo Alto, CA 

94304 , USA 

'^Institute of Theoretical Astrophysics, University of Oslo, P.O. Box 1029 Blindern, N-0315 

Oslo, Norway 

^National Solar Observatory, PO Box 62, sunspot, NM 88349, USA 
^ Centre for Mathematical Plasma Astrophysics, Mathematics Department, KU Leuven, 
Celestijnenlaan 200B bus 2400, B-3001 Leuven, Belgium 
^ The University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
^Current address: STEL, Nagoya University, Aichi 464-8601, Japan 

Patrick.antolinOnao.ac.jp 

ABSTRACT 


Transverse magnetohydrodynamic (MHD) waves are ubiqnitons in the solar 
atmosphere and may be responsible for generating the Sun’s million-degree outer 
atmosphere. However, direct evidence of the dissipation process and heating from 
these waves remains elusive. Through advanced numerical simulations combined 
with appropriate forward modeling of a prominence flux tube, we provide the ob¬ 
servational signatures of transverse MHD waves in prominence plasmas. We show 
that these signatures are characterized by thread-like substructure, strong trans¬ 
verse dynamical coherence, an out-of-phase difference between plane-of-the-sky 
motions and LOS velocities, and enhanced line broadening and heating around 
most of the flux tube. A complex combination between resonant absorption and 
Kelvin-Helmholtz instabilities (KHI) takes place in which the KHI extracts the 
energy from the resonant layer and dissipates it through vortices and current 
sheets, which rapidly degenerate into turbulence. An inward enlargement of the 
boundary is produced in which the turbulent flows conserve the characteristic 
dynamics from the resonance, therefore guaranteeing detectability of the reso¬ 
nance imprints. We show that the features described in the accompanying paper 
( Okamoto et al.||2015 ) through coordinated Hinode and IRIS observations match 
well the numerical results. 
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Subject headings: magnetohydrodynamics (MHD) — Sun: activity — Sun: corona 
— Sun: filaments, prominences 


Introduction 


Observational reports of periodic transverse displacements of chromospheric and coronal 
structures have increased in the last decade, especially due to the advance in instrumentation 
allowing higher sensitivity and spatial, temporal and spectral resolution. Initially reported 
through EUV imaging as large disturbances of coronal loops in the aftermath of flares with 
TRACE (Nakariakov et aL]|1999), they are now observed more frequently with SDO (McIn¬ 


tosh et ah 2011; Threlfall et al. 2013; Thurgood et al. 2014). Observations with the Solar 


Optical Telescope (SOT, [Tsuneta et al.|2008 ) on board Hinode (Kosugi et al.|2007), together 
with other ground-based facilities providing high spatial resolution imaging, have established 


that such low amplitude periodic disturbances are ubiquitous 

in the chromosphere (De Pon- 

tieu et al. 

2007b|a; 

Kuridze et al. 

2012; 

Morton et al. 2012 

) and 

in the corona, through 

observations of prominences or coronal rain (Okamoto et al. 2007 

; Terradas et al. 2008b; 

Lin et al. : 

2009; Lin 2011; Antolin & Verwichte 20111 |Arregui et al. 2012; Hillier et al. 2013; 

Schmieder et al. 2013). Periodic sj 

rectral disturbances in Doppler velocities have also been 
h corresponding swaying in intensity of chromospheric or 

reported, frequently associated wit! 

coronal structures ( 

Tomczyk et al. 

2007; 

Erdelyi & Taroyan 

2008; 

Van Doorsselaere et al. 

2008; Tomczyk & McIntosh 2009; ! 

Lin et al. 2009; Tian et al. 

2012; 

Schmieder et al. 2013). 


The ubiquity of such motions and also of torsional components has recently been further es 
tablished in spicules, by the Interface Region Imaging Spectrograph {IRIS, |De Pontieu et al. 

mibff 


Rouppe van der Voort et al. 2015; Skogsrud et al. 2014) 


Such periodic transverse disturbances usually have low amplitudes (a few km 


.-I'l 


periods in the order of a few to tens of minutes and usually exhibit damping signatures. 
General consensus exists on the interpretation of such atmospheric disturbances as Alfvenic 
waves, specifically transverse MHD waves (e.g., Oliver j20'09; Mackay et al.|2010; Arregui et al. 


2012; Mathioudakis et al. 2013, and references therein), also known as kink waves or kink 


surface Alfven waves (Goossens et al. 2002, 2012). The importance of these Alfvenic waves 


is set partly on the close connection between their characteristics and the medium in which 


they propagate, thereby serving as seismological probes of the magnetic field (Nakariakov 


& Ofman 2001), the loop radial structure (Aschwanden et al. 2003), the vertical density 


scale height (Andries et al. 2005) and Alfven transit times (Arregui et al. 2007). Alfvenic 


waves are further important due to their ability to carry significant amounts of energy over 
long distances, thereby potentially playing an important role in coronal heating and solar 
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wind acceleration (Alfven 1947; Uchida & Kaburaki 1974; Antolin & Shibata 2010 


van 


Ballegooijen et al. 2011; Matsumoto & Suzuki 2014). These waves are often observed to 


damp on short spatial and temporal scales ( Aschwanden et al.||1999 i Pascoe et al.j2010 , 2011). 
However, it is still unclear whether the damping usually leads to signihcant dissipation, and 


thus heating (see Parnell & De Moortel 2012, for a contemporary review). 


The dissipation process of transverse MHD waves in an inhomogeneous medium has 

Ionson||1978 ; Hollweg et al.||1990 


been a subject of active research over several decades (e.g. 


Sakurai et al. 1991; Goossens et al. 1992). Theoretically it is expected that the transverse. 


bulk motion is subject to damping by resonant absorption in an inhomogeneous flux tube. 


such as a prominence thread (Arregui et al. 2008; Arregui & Ballester 2011; Soler et al. 


2012). An analogous mechanism known as field line resonance exists in the magnetosphere 
(Southwood &: Hughes 1983). Resonant absorption occurs at the tube’s boundary layer. 


where the phase speed of the trapped (anisotropic) kink wave matches that of the external 
(isotropic) Alfven waves. The kink wave (with azimuthal wave number m = 1), whose initial 
energy is mostly transverse, mode converts into azimuthal waves locally resembling torsional 
Alfven waves, which are finely localized around the tube’s boundary layer ( Verth et al.||2010 ; 
Arregui et al. 2011). Azimuthal motions are thus amplified and introduce velocity shear. 


making them prone to be unstable to the Kelvin-Helmholtz instability (KHI) (Heyvaerts & 


Priest 

1983; 

Soler et al. 

2010 

). The KHI deforms the boundary layer and leads to enhanced 


dissipation of the wave energy into heat in thin, turbulent current sheets (Uchimoto et al. 


1991[|Karpen et al.|1993[|Poedts et al.|1997[ [Ziegler fc Ulmschneider|1997t [Otman et al.|1998 


Terradas et al. 2008a). Apart from chromospheric and coronal heating, the importance of 


this instability has been stressed in magnetic reconnection studies (Lapenta & Knoll 2003) 


and in the mixing of plasma between two magnetically confined plasma environments, such 


as closed and open magnetic field lines (relevant for solar wind studies Andries & Goossens 
2001) or across the magnetospheric boundary ( [Fujimoto fc Terasawa 1994; Fujimoto et al. 
200 ^. 


Recent simulations by Antolin et al. (2014) show that even small amplitude oscillations 


can lead to the KHI, whose vortices combined with line-of-sight (LOS) effects result in 
strand-like structures within coronal loops, but at scales that are difficult to resolve with 
current coronal instruments. This study opens the question of whether an important part 


of the observed strand-like structure in loops (< 

Dfman & Wang 2008; Reale 

: 2010; Antolin 

& Verwichte 2011; Antolin & Rouppe van der Voort 2012[ [Brooks et al. 2013; Peter et al. 

2013; 

Scullion et al.||2014) or thread-like structure in prominences (Lin||2011; 

Okamoto et al. 

2007; 

Terradas et al.||2008b; Vial & Engvold||201f 

)) could be linked to transverse MHD waves 


through their ubiquity and the generation of the KHI. 
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In this series of two papers, unique coordinated observations with i/mode/SOT of Hin- 
ode and IRIS are used to provide evidence for a telltale sign of resonant absorption and 
associated KHI and heating in an active region solar prominence ( Qkamoto et al.|2015 i here¬ 
after Paper 1). In the present paper details of the three-dimensional MHD simulations and 
forward modeling first introduced in Paper 1 are given, and we further provide a complete 
description of the observational signatures expected for transverse MHD waves. As discussed 
in this paper most of the obtained signatures in the numerical model are also expected in 
other solar atmospheric structures such as coronal loops (Antolin et al. 2014) and spicules, 
and may be considered general to a certain extent. 


The paper is organized as follows. In section the details of the numerical model are 
provided. The results of the forward modeling of the numerical results are given in section 
Results are discussed in section |4] and conclusions in section [5l 


2. Numerical Model 


To interpret the observations in Paper 1 we performed 3D MHD simulations of a promi¬ 
nence fiux tube oscillating with a fundamental kink mode. The numerical model is similar 
to that of the previous work Antolin et al. (2014) for the coronal case, and we adapt it 
to the prominence case. We considered a fiux tube with physical quantities usually found 
in solar prominences (Tandberg-Hanssen |1995| [Vial &: Engvold 2015). We take a number 
density and temperature contrast of 10 and 1/100, respectively, with respect to the external 
values, which we set to 10® cm“^ (10^^ m“^) and 10® K. The external plasma is assumed 
to be fully ionized with a beta value of 0.01, leading to an external magnetic field value of 
18.6 G (1.86 mT) and a slightly larger value within the fiux tube in order to keep a total 
pressure balance. The external and internal thermodynamic values are linked smoothly as 
in Antolin et al. ( 2014| ), defining a specific width I for the boundary layer (see Fig. [^. Here 
we adopt a width of l/R ~ 0.4, where R denotes the radius of the tube. The dimensions 
of the flux tube are to some extent constrained by the observations. AIA reveals that the 
length of the prominence is larger than 150,000 km and exists over a height range of roughly 
30,000 km. However, a strong LOS projection effect is expected. Also, dynamic coherence 
is only detected over a few 1,000 km width in the vertical direction (fig. 2 in Paper 1), 
suggesting that the prominence is likely to be composed of several dynamically independent 
fiux tubes. We thus choose to take the length of the tube as 200 R and we set i? = 1,000 km. 


At time t = 0 the fiux tube is given a perturbation in velocity mimicking the fundamental 
kink mode (with longitudinal wavenumber k ~ 0.015/R) as described in Antolin et al. (2014). 
The amplitude of the pulse is Uq ~ 8 km s“^, a value in the low end of those being detected 


Antolin et al. (2014 
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Fig. 1.— Sketch of the prominence flux tube subject to a transverse MHD wave in our 
numerical model. The axis of the flux tube is oriented along the 2 ; axis. The transverse 
perturbation (and the subsequent oscillation) is directed along the a:-axis. The LOS for the 
forward modeling is in the (x, y) plane. The 0° and 90° LOS rays are directed along the 
positive y — axis and x-axis, respectively, as shown in the figure. Variables p, T, p and B 
denote, respectively, the density, temperature, pressure and magnetic held within (subindex 
i) or outside (subindex e) the flux tube. R and I denote, respectively, the radius of the flux 
tube and the width of the inhomogeneous boundary layer. 


with SOT and IRIS. The kink phase speed is Ck ~ 776 km s and therefore the perturbation 
is largely linear. 


The numerical simulations were performed with the CIP-MOCCT scheme (Kudoh et al. 


1999), which solves the 3D MHD equations including resistivity and viscosity. Several simu¬ 


lations were performed varying the values of resistivity and viscosity, and the number of grid 
points. The base run consists of a numerical box with 512 x 256 x 50 points in the x, y and 
.2 directions, respectively, where x and y denote the parallel and perpendicular directions 
of oscillation, respectively, and z is the longitudinal direction. Higher resolution runs with 
double number of grid points in every direction were also performed. The values of resistivity 
and viscosity are set to small constant values everywhere, ensuring large enough Lundquist 
(and viscous) Reynolds numbers of 10^ — 10^ in the base run, and around 10® — 10^ in the 
large resolution run. While a considerable increase of fine scale structure is present in the 
high resolution runs, the main results are left unchanged. This is expected since the effects 


from resonant absorption stabilise for high Reynolds numbers (Kappraff &: Tataronis 1977 
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Poedts & Kerner 1991 

; Ofman & Davila 19951 |Tirry & Goossens 1996 

; Van Doorsselaere 

et al. 

2004; 

Terradas et al. 

2006 

). The results presented here are from the base run unless 


explicitly stated. 


Thanks to the symmetric properties of the kink mode only half the plane in y and z is 
modeled, and we set symmetric boundary conditions in all boundary planes except for the 
X boundary planes, where periodic boundary conditions are imposed. In order to minimize 
the influence from side boundary conditions the spatial grids in x and y are non-uniform, 
with exponentially increasing values for distances beyond the maximum displacement. The 
maximum distance in x and y from the center is ~ 16 R. The spatial resolution at the loop’s 
location is 0.0156 R, and a tenth of this value for the large resolution run. 


2.1. Resonant absorption in the numerical simulation 


Following the initial kink the tube oscillates with a period of 525 s, close to the funda¬ 
mental period of 2L/ck- The maximum displacement of the tube is 710 km and the damping 
time of the oscillation is roughly 840 s, a value 16% smaller than that predicted by resonant 
absorption theory in the linear regime. As the tube is displaced transversely by the kink 
mode resonant absorption sets in immediately, as expected from theory ( Poedts et al.||1990 ). 
This can be seen in movie 1, as black and white rings in the z-vorticity panel, occurring 
around the boundary layer. In the boundary layer, where the Alfven speed is close to the 
kink speed, the nature of the kink quasi-mode is modified and becomes similar to that of a 
torsional Alfven wave. Energy is transferred from transverse motions to azimuthal motions 
around the resonant layer, seen as bright or dark rings in the z-vorticity panel. The process 
of resonant absorption has been extensively studied in the past and has been shown to be 


(Goossens et al. 

2002, 

2006; 

Van 


Doorsselaere et aL]|2004 ; Goossens et al.||2011 ; Pascoe et al.||2010 ). It has been further shown 
to be a robust mechanism, with little dependence on the degree of density inhomogeneity 


within the tube (De Moortel & Nakariakov 2012; Pascoe et al. 2011; Terradas et al. 2008c). 


2.2. Kelvin-Helmholtz feedback from resonant absorption 

Apart from resonant absorption, the KHI sets in rapidly, just one period after the initial 
kink, as seen in movie 1. The instability is generated by the shear fiows at the edges of the 
loop, produced by the kink mode (from the dipole-like azimuthal flow outside the tube 
and the transverse displacement of the tube). Such flows have been shown to generate the 
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Fig. 2.— Snapshot of the cross-section of several physical quantities at the apex of the 
prominence thread at a time in the simulation of t = 1162 s, corresponding to roughly 2 
oscillation periods. From top to bottom the panels in both columns correspond to number 
density, temperature (logarithm), flow velocity fleld, and the longitudinal (z-) components 
of the vorticity and current density. The left and right panels correspond, respectively, to 
the base and high resolution runs. An animation of this figure for the high resolution run is 
available online (movie 1). See also movie B in Paper 1 for the corresponding animation of 
the base run. 


instability in photospheric and coronal conditions (Karpen et ah 1993; Ofman et al. 

1994; 

Poedts et al.||1997 

Ziegler & Ulmschneider||1997; Terradas et al.||2008a), for small amp] 

itude 

n the 

kink modes (Antolin et al. 2014) matching the currently observed transverse motion i 

solar atmosphere ( 

Okamoto et al. 2007; Tomczyk et al. 2007; Erdelyi & Taroyan 2008 

i Van 

Doorsselaere et al.|2008i Terradas et al.|2008bi Lin et al. |2009| |Lin|2011i McIntosh et al. 

2011; Tian et al. ! 

2012; Antolin & Verwichte||2011: |Arregui et al.||2012l |Hillier et al. 

2013; 

Schmieder et al. 2013). Resonant absorption helps the onset of the instability by increasing 


the amplitude of the azimuthal flow in the boundary layer and therefore that of the shear flow. 
Here we show that the instability can easily be obtained in prominences. The fast growth of 






































































































the unstable modes is guaranteed by the large density contrast and the small longitudinal 


wavenumber, as shown in Antolin et al. (2014). In the present case, application of formula 2 


in the previous paper predicts the most unstable mode with azimuthal wavenumber m = 3. 
In the base and high resolution run 3 and 5 large eddies, respectively, are initially obtained, 
close to the theoretical results from linear analysis. 

The instability produces vortices (eddies in 2D cross-sections) that develop in the plane 
perpendicular to the axis of the flux tube, as shown in Fig. Such coherent structure 
rapidly degenerates in turbulence as small scales are produced from the continuous shear 
flows. The vortices’ lifetime is about 1 period, and they are continuously generated and 
destroyed during the entire simulation. It is important to notice that the vortices do not 
only affect the boundary of the flux tube but also penetrate deeply towards the core. This 
can be seen clearly in movie 1 and Fig. As a result, the ensuing turbulence enlarges the 
boundary layer and rapidly drifts inwards towards the center of the tube. Strong deformation 
of the tube’s cross section is clearly visible within one period from the onset of the KHI: 
the initially dense cool core is significantly reduced to roughly half its original size, and the 
boundary layer has enlarged to roughly the same size of the core (l/R ~ 1). The enlarged 
boundary layer has a temperature around 10^ K, meaning that most of the tube has been 
subject to heating. The nature of this heating is not only viscous but also comes from 


ohmic dissipation. Indeed, we obtain current sheets at the vortices (Ofman 2009; Antolin 


et al.||2014 ), therefore spanning along most of the prominence flux tube, due to the relatively 
strong shear in the magnetic field. 

Component magnetic reconnection can occur at such current sheets, as has previously 


been shown (Lapenta & Knoll 2003). The heating that ensues from viscous and ohmic 


dissipation can be strong and impulsive at the start of the KHI exhibiting various short 
lived bursts, but rapidly sets in to a constant and continuous tendency. The KHI effectively 
extracts the kinetic and magnetic energy in the boundary layer deposited by the resonant 
absorption mechanism, and converts it into heat through viscous and ohmic dissipation from 
the turbulence and current sheet generation. The regions around the resonant layer affected 
by the KHI heat up rapidly. This means that most of the tube is heated from chromospheric 
to transition region temperatures in the time scale of one period. In the high resolution 
run the deformation and inward enlargement of the boundary layer is more severe. This is 
because of the enhanced turbulence due to the significant increase of fine-scale structure. 
The number of vortices and current sheets increases leading to wider and faster heating from 
a more efficient energy cascade towards smaller scales. The timescale of the heating of most 
of the flux tube to transition region temperatures is decreased by a few minutes compared 
to the base run. Only the inner core is more slowly and gradually heated. It is important 
to notice that the energy transfer from the resonant layer towards the inner parts of the 
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flux tube in our model does not depend on the perpendicular thermal conduction, which 
is expected to be very small even in the case of prominences. This transfer is produced 
by the constant generation of small-scale current sheets and vortices in the turbulent flows 
that occur in the perpendicular magnetic held direction (but which never cross the magnetic 
held). 


2.3. Resistivity in the numerical model 


The Lundquist number in the base run is on the order of 10^ — 10^, and 10® — 10^ in 
the high resolution run. On the other hand, the resistivity in the solar corona is expected 
to be many orders of magnitude less than the numerical resistivity, indicating much higher 
Lundquist numbers in the solar corona than those considered here. However, the development 
of the KHI and phase mixing leads to the rapid formation of gradients on smaller and smaller 
scales thereby lessening the impact of the resistivity issue. In addition the prominence 
consists of chromospheric material in which neutral particles naturally occur. As a result, 
the resistivity in chromospheric plasmas can be on the order of that in numerical models. 


especially in the presence of ambipolar diffusion (Martinez-Sykora et al. 2012). 


2.4. Limitations of the numerical model 


In this numerical model gravity is not included, implying no density stratification along 
the flux tube. Thermal conduction and radiation are not included either. On the geometrical 
side, the flux tube is straight, hence effects from curvature, expansion or coherent twist are 
also neglected. Finally, longitudinal flows are not taken into account either. While each of 
these factors is present in the observed prominence to a certain degree, we believe that they 
do not contribute substantially to the physical mechanisms discussed here. Each of these 
mechanisms has first order contributions in the longitudinal direction and are only second 
order in the transverse direction, which is where the important dynamics for the present 
investigation take place. Furthermore the timescales of resonant absorption or the KHI are 
significantly faster than those of thermal conduction or radiation. Prominences such as the 
one observed here are mostly horizontal with respect to the solar surface, implying a minor 
role for stratification, expansion or curvature. Flows can play an important role on damping 
or enhancement of transverse oscillations (Gruszecki et al.||2008; Terradas et al.||2010), which 


can limit the efficiency of resonant absorption and the onset of KHI (Soler et al. 2012). 


However, we consider this effect not significant due to the very rapid timescales of these 
mechanisms. 
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Another limiting factor of our model is the symmetry, especially in the r/-direction, 
based on the spatial properties of the ideal kink mode. This assumption implies that no 
changes in polarisation nor coupling between azimuthal wave modes take place, and that 
therefore most of the energy will remain in the m = 1 kink mode polarised along the 
X—direction. Density inhomogeneities, like those produced from the KHI in our model, can 
lead to mode coupling. Such effects, which we consider secondary, may lead to detectable 
Doppler asymmetries between diametrically opposite regions in the flux tube, and should be 
studied in more detail in the future. 


The symmetry of our model implies the absence of twist, which may probably be the 
most limiting factor, but only on the KHI. A small amount of twist may inhibit the instability 
in the linear regime (Soler et al. 2010). However, the presence and degree of twist in a 
prominence is a subject of debate. Also, no numerical parameter study exists investigating 
the full effect of twist in a more realistic model. On the other hand, a transverse MHD wave 
still produces multiple current sheets in a flux tube with full coherent twist Ofman (2009). 
In our simulations the KHI produces partial twist of the flux tube, which slows down the 
instability itself but does not suppress it. 


3. Forward modeling 

3.1. Conventions on LOS angle and spatial resolution 

For the forward modeling we define the LOS angle such that 0° is parallel to the positive 
y-axis, and 90° is parallel to the positive x-axis, which is the axis of oscillation (see Fig. [^. 

For correct comparison with observations with a given instrument with resolving power 
of X we degrade the original spatial resolution of the numerical model by first convolving 
the image of interest with a Gaussian with FWHM of X, simulating a PSF of roughly the 
same size. For each pixel we then sum over a region of X in size centered on the pixel. 

When comparing with observations it is important to keep in mind that the spatial 
resolution of SOT (0.2") is approximately twice as good as that of IRIS (0.33" — 0.4"). 
Accordingly, the imaging and spectroscopic results presented here have been convolved with 
different Gaussian functions whose FWHM ratio is in accordance to the PSF ratio between 
SOT and IRIS. This is especially important when predicting the emission in slits that are 
relatively far away from the observed threads, as in the present case. 
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3.2. Optically thick approximation 


Radiative transfer calculations have been performed for synthesis of the Mg 11 h&k and 
Ca II H&K resonance lines in non-LTE (see Paper 1). We used the two-dimensional version 


of the RH code (Uitenbroek 2001), based on the method described in Rybicki & Hummer 


(1992). Results with RH indicate that only rays crossing the prominence core are optically 


thick and that only the boundary layer contributes to the emerging intensity. Furthermore, 
the intensity in the boundary layer is optically thin, and therefore mostly temperature depen¬ 
dent. Since RH is computationally intense, we have adopted an approximation to optically 
thick emission with the FoMo code (https://wiki.esat.kuleuven.be/FoMo, [Antolin h 
Doorsselaere||20 13 ), which allows to perform forward modeling of time dependent results for 
the Mg II k line for the entire numerical box. This approximation is based on the RH results 
described in Paper 1 and consists solely of limiting the integration path of the emission line 
flux (G\{T,n)nl function) along a given LOS (T and rig denote, respectively, the tempera¬ 
ture and electron number density at a given position along the LOS), as follows. At a given 
time t during the simulation, the specific intensity at wavelength A along a given LOS is 
proportional to: 


I such that 


rOT 

X,ray 


(X 



Gx{T, ne)nldl', 


( 1 ) 


pl pLmax 

/ Gx{T{t),nXt))nXtfdl' = OTx / Gx{T{Q),nM>e{^fdl' ( 2 ) 

Jo Jo 

where Lmax is the maximum path length along the LOS and OT is a parameter between 0 
and 1 indicating the degree of optical thickness. OT= 1 is the optically thin resullj^ The 
left hand side is the integral up to a distance I along the LOS of the emission line flux (and 
therefore part of the total intensity). The limit of the integral, /, is physically similar to the 
distance along the LOS where the optical thickness at the wavelength A becomes unity. The 
right hand side of the equation is the optically thin intensity at time t = 0 (disregarding the 
abundance of the element) multiplied by the factor OT. Fixing the parameter OT defines 
I, and therefore the unknown opacity of the material is determined. We therefore assume 
that the optical thickness is roughly inversely proportional to the optically thin intensity of 
the material. As a proportionality factor we assume that the maximum achievable intensity 
along such LOS rays is roughly the optically thin intensity at time t = 0 (right hand side 
with OT= 1). This is reasonable since the model stays mostly in the linear regime (no major 


^this is the case for t = 0, and is a sufficiently good approximation for t > 0 
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longitudinal perturbation) which leads to little mass displacement along the flux tube axis. 
For instance, with a value of 0T= 0.5 for a LOS crossing the center of the flux tube only 
rays starting from the center contribute to the emerging intensity at time t = 0. 

In optically thin conditions the emission line flux in the Mg II k line obtained from 
CHIANTI has a peak at log(r) ~ 4.2, leading to a ring shaped structure similar to that 
of the source function obtained from RH (see Fig. 8 in Paper 1). This optically thick 
approximation therefore recovers the essential aspects of the radiative transfer calculations 
and is justified by the mainly temperature dependent source function in Mg II k and Ca II H. 


3.3. Thread-like structure, out-of-phase behaviour, edge broadening and 

heating 


In the top left panel of Fig. [^we show the optically thin case of the Mg II k intensity 
for a theoretical slit located perpendicular to the flux tube, at the apex, and for a LOS 
angle of 45°. Here we further assume that the slit has a spatial resolution equal to that of 
the numerical resolution. As for the coronal case (Antolin et al. 2014), thread-like structure 
rapidly appears after just one period of oscillation in the prominence case and permeates 
the entire flux tube. The threads correspond to the KHI vortices, which, as explained in 
section |2.2| initially occur at the surface of the flux tube boundary and rapidly degenerate 
in turbulent-like flows drifting towards the core of the prominence. 


From the synthetic emission profiles obtained with FoMo we calculate the moments by 
single Gaussian fits to the line profiles (double Gaussian fits are not necessary due to the 
relatively low velocity amplitudes). The middle and lower left panels of Fig. [^correspond, 
respectively, to the LOS velocity and line broadening obtained in this way from the previ¬ 
ously mentioned theoretical slit. In the Doppler velocity panel, apart from the initial kink 
displacing the flux tube in the POS, the onset of resonant absorption can be seen at the 
edges right from the beginning as enhanced Doppler signals, which become rapidly shifted 
in time due to phase mixing (more precisely, due to the change of the Alfven speed across 
the boundary layer). This phase shift of the azimuthal flows at the flux tube boundary is 
ultimately responsible for the out-of-phase behaviour with the POS motion, as explained in 
Paper I, which we further elaborate in the following sections. Resonant absorption is also 
visible in the line broadening panel, where clear enhancement is highly localized at the flux 
tube boundaries, where the azimuthal flows are directed along the LOS. 


As the KHI is triggered, at this high spatial resolution traces of the KHI vortices can be 
distinguished in the x—t diagrams of the Doppler velocity and line broadening like thread-like 
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structure similar to that in the intensity panel. This reflects the fact that the vortices, while 
embedded in the large-scale flow generated by the kink wave, have their own characteristic 
small-scale vortex flows. The KHI extracts the energy from the resonant layer and dumps it 
into heat through ohmic and viscous dissipation at the generated current sheets and vortices. 
A result of this mechanism is the generation of turbulence as the energy of the large eddies 
cascades down to lower eddies. This is clearly shown in movie 1 (see also movie 2 in Paper 1), 
where the higher resolution allows the generation of a higher number of small-scale current 
sheets and vortices, leading to a higher rate of dissipation. This results in an enlargement of 
the flux tube boundary, where significant heating occurs, and a corresponding enhancement 
and enlargement of the line broadening around the flux tube boundary, as shown in the 
lower left panel of Fig. Importantly, this enlargement is absent in the Mg II k intensity 
panel but is present in the Doppler velocity maps. This effect will further be analysed in 
section 1331 

The heating can also be clearly observed in synthetic Si IV images, as shown in the right 
column panels of Fig. This column shows the same quantities as in the left column but 
for Si IV 1402.77 A, which has a formation temperature around 80,000 K. As the flux tube 
with its strand-like structure fades out in Mg II (top left panel), similar structure fades in in 
Si IV (top right panel). The spectral signatures shown in the middle and lower right panels 
(Doppler velocity and line broadening, respectively) have essentially the same features as 
for Mg II. However, it is worth noting that the fine-scale structure in intensity, Doppler and 
line width is more clearly seen in the TR line at this high resolution, since the temperature 
formation of this spectral line captures better the heating events (current sheets and vortices) 
occurring around the boundary. 


3.4. Effect of optical thickness 

In the simple geometry considered of a single flux tube we show in Fig. 7 of Paper 1 
that chromospheric emission in Mg II h&k and Ca II H&K comes from the most dynamically 
important (ring-shaped) section of the flux tube, where resonant absorption and dynamic 
instabilities such as KHI take place. The ring shape of the source function in Mg II k and 
Ca II H are a temperature effect from the optically thin layer of material surrounding the 
prominence core. These results justify the approximation to the Mg II h&k intensity based 
on optically thin calculations with the FoMo code explained in section m The result of 
this approximation is shown in Fig. where the emerging intensity in Mg II k for several 
values of the OT parameter is plotted with respect to time for a theoretical slit located at 
the flux tube apex, perpendicular to the axis and for a LOS of 45° and a spatial resolution 
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LOS velocity for 45° LOS 



Si IV 1402.77 optically thin intensity for 45° LOS 



Line broadening for 45° LOS 
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Time [min] 


Line broadening for 45° LOS 
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Fig. 3.— From top to bottom we show time slices of the optically thin intensity, LOS velocity 
and line broadening for the Mg II k line (left panels) and the Si IV 1402.77 A line (right 
panels) that a theoretical slit would capture being perpendicular to the prominence axis, 
located at the apex and making a LOS angle of 45° (see section 3.1). We further assume 
that the slit spatial resolution is equal to the numerical resolution. LOS velocity and line 
broadening are calculated with single Gaussian htting to the line profiles. 


of 0.1 R. The decrease in intensity with decreasing value of the OT parameter is due to 
the lower column density contributing to the emerging LOS ray. For instance, the case of 
OT = 0.5 corresponds to an optically thick case in which the radiation along the center of 
the tube crosses roughly half the tube only. Several values of the OT parameter have been 
tried, from 0.1 to 0.7. We can see that thread-like structure is always obtained (except the 
case OT= 0.1, which is too dark relative to the optically thin case), although the number of 
threads increases for decreasing optical thickness, as expected. 

In Fig. [^we show the LOS velocities that result for the optically thick case OT= 0.5 
of Mg II k. The figure corresponds to the time slice from a theoretical slit located at the 
center of the tube, perpendicular to its axis, for different LOS angles. The effect of optical 
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thickness on LOS velocity is only significant when observing at shallow LOS angles. In other 
words, optical thickness in general does not affect the measured LOS velocity. Indeed, due 
to the symmetry of the kink mode, only at low angles (< 15°) the Doppler signals from 
opposite regions with respect to the loop’s center plane (x-axis) tend to cancel each other. 
If the loop’s core is opaque in a given wavelength only one side of the loop’s surface will 
contribute, avoiding cancellation. In Fig. Doppler signals can be seen even for a 0° angle. 


Intensity Mg II 2796.35 Opt. thick approx. Resolution = 0.1 R 45° LOS plane 
Optically thin case 



0 10 20 30 40 50 

Time [min] 


Fig. 4.— Effect of the optically thick approximation explained in section on a time slice 
of the emerging Mg II k intensity for a slit perpendicular to the prominence axis, located at 
the apex and making a LOS angle of 45°. From top to bottom we show increasing degrees 
of the optically thick parameter OT, from purely optically thin (OT= 1) to highly optically 
thick (OT— 0.1). We show the effect for a spatial resolution of 0.1 R. 


3.5. Effect of LOS angle 


In Fig. I^we show the effect of LOS angle on the emerging intensity in Mg II k for the 
case of optically thin emission and a spatial resolution of 0.1 R. As seen in the figure, the 
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Doppler velocity Mg II 2796.35 OT = 0.5 Resolution = 0.0 R 
Variation along 0° LOS piane 



Variation aiong 15° LOS piane 



Variation aiong 30° LOS piane 



Variation aiong 45° LOS piane 



Variation aiong 60° LOS piane 



Variation aiong 75° LOS piane 



Variation aiong 90° LOS piane 



Time [min] 


[km/s] 


5 
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Fig. 5.— Time slice of the Doppler velocity for Mg II k along a slit perpendicular to 
the prominence axis, located at the apex. The parameter OT is set to 0.5, corresponding 
to a moderately optically thick case. The spatial resolution is the same as the numerical 
resolution. Each panel corresponds to a different LOS angle. 


thread-like structure appears at any angle. Due to LOS effects, apparent crossing of threads 
is observed, especially for the optically thin case. This corresponds in the model to either 
twisting within vortices or to the combination of several vortices along the LOS. As in the 
coronal case the number of threads and the distribution of widths is highly dependent on the 



17 


Intensity Mg II 2796.35 Optically thin case Resolution = 0.1 R 
Variation along 0° LOS piane 



Variation aiong 15° LOS piane 



Variation aiong 30° LOS piane 



Variation aiong 45° LOS piane 



Variation aiong 60° LOS piane 



Variation aiong 75° LOS piane 



Variation aiong 90° LOS piane 
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Fig. 6.— Time slice of the optically thin intensity for Mg II k along a slit perpendicular 
to the prominence axis, located at the apex. The spatial resolution is 0.1 R. Each panel 
corresponds to a different LOS angle. 


initial kink amplitude and the width of the boundary layer. Similarly, since vortices move 
coherently (they are embedded in the large scale azimuthal flow) threads detected at coarse 
spatial resolution actually correspond to groups of threads, as can be seen comparing Figs. 
and|^ The effect of spatial resolution will further be discussed in section |4T 
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The preferential occurrence of resonant absorption and KHI vortices along the direction 
of oscillation (with azimuthal flows preferentially directed along the x-axis) generates an 
asymmetry in the flux tube cross-sections for all physical quantities, making these more 
elliptic. Indeed, the generation of current sheets and vortices where most of the heating 
takes place is concentrated in lobes around the prominence core, with maxima at ±90° 
from the a:-axis and minima along the a:-axis (clearly seen in movie 1). These lobes are 


expected from the spatial distribution of the m = 1 mode (as shown also by Pascoe et al. 


2011, ; see their Fig. 4). A corresponding decrease of density and increase of temperature is 
obtained along these lobes, which broadens inwardly the flux tube boundaries. The emission 
line flux in Mg II k is thus diminished, partially due to the density squared dependence of 
the intensity, which is expected under the optically thin conditions of the boundary layer 
(contributes in generating sharp threads), but especially to the increase in temperature in the 
same layer, which brings the plasma gradually out of the temperature formation of Mg II k. 
Accordingly, this effect leads to a thinning of the flux tube in time for most LOS angles, 
as shown in Figs. [^and[^ Also, a difference in width of the flux tube is readily seen when 
observing at different LOS angles, as shown in Fig. On the other hand, the thinning 
is not observed in synthetic Si IV intensity images, as shown by Figs. and since the 
temperature of the heated plasma enters the range of temperature formation of this spectral 


line. This LOS angle effect has also been reported in the coronal case (Antolin et al. 2014). 


Through cross-correlation analysis we And that the lifetime of a KHI vortex is roughly 


one period. However, as in the coronal case (Antolin et al. 2014) it is diflicult to estimate 


precisely in the model the lifetime of threads, due to the continuous generation and destruc¬ 
tion of vortices and the projection effects. A new vortex can appear at the location of a 
previously destroyed vortex, thus producing apparently longer lived thread-like structure. 
This can be checked in Fig. and movie 1 (especially clear in the high resolution run). Fur¬ 
thermore, we can expect such vortices to be longer lived in the real corona since turbulence 
is only independent on viscosity for larger viscous Reynolds numbers. 


The effect of the LOS angle on the LOS velocity is shown in Fig. where the optically 
thick case of OT= 0.5 is shown for Mg II k. The time slice clearly shows that the red- 
blue patterns at the edges of the tube are shifted towards 90° out of phase with respect to 
the Doppler shifts at the center, and 180° out of phase with respect to the transverse tube 
displacement in the POS. This results in characteristic arrow (>) shaped structures in such 
time slices, which can be observed for basically any LOS angle and any optical thickness 
(even at shallow angles as long as the line opacity is significant, as explained before). This 
characteristic out-of-phase behavior from resonant absorption can be clearly seen after just 
half an oscillation period (before the onset of the KHI), and remains regardless of the presence 
or absence of the instability. However, the detectability of this out-of-phase behavior and 
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the instability is significantly increased because of the spatial broadening of the boundary 
layer that the instability produces. 


In Fig. [^we show the effect of LOS angle on the measured line broadening of the Mg II k 
line. As for the LOS velocity, the maximum values are found around the boundary layer, 
highly localized at first in the resonant layer before the occurrence of the KHI, and more 
widespread after occurrence. Part of this line broadening is not due to resonant absorption 
but to the fact that at the edges the LOS crosses more material with different velocities 
(especially in optically thin conditions), which can clearly be seen at t = 0 for a 0° LOS 
(angle for which the azimuthal component is largely absent at the edges). This effect has 
also been described by Doorsselaere & Nakariakov (2008). As seen in the figure, fine broad¬ 
ening becomes more detectable (occupies an apparently larger area of the flux tube) when 
observing along the axis of oscillation. This is expected since such LOS rays cross the lobes 
of oscillation of the m = 1 mode. The maximum values, however, are detected at shallow 
angles, just after KHI occurrence. This is because LOS rays at such angles cross vortex 
flows with different vorticity occurring next to each other. The flow within a given vortex 
introduces strong velocity shear, especially when embedded in the azimuthal flow amplified 
by the resonance. Further small vortices are then generated on top of the initial vortex, 
with different small-scale flow direction (different vorticity) but still embedded in the larger 
azimuthal flow. This can be clearly seen in the vorticity panels of Fig. and in movie I 
(see also movie 2 in Paper I). These results could potentially explain observations of line 
broadening enhancement at loop edges or in inter-loop planes (Doschek et al. 2007). 


4. Discussion 
4.1. Intensity variability 

As discussed in Paper 1 two kinds of intensity variability can be noticed in the obser¬ 
vations of prominence threads. The first kind is linked to flows or waves in the threads and 
occurs on short timescales of one oscillation period or less. The second kind is linked to the 
dimming in Ca II H intensities and corresponding intensity increase in Si IV and is therefore 
linked to heating. The latter is more gradual and occurs in timescales of I — 2 periods. 

As in the observations, the intensity variability obtained in our simulations has also 
two components: a short timescale component linked to dynamics and a slightly longer 
component linked to large scale heating. These two components can be seen in Fig. (also 
Figs. 1^ and 1^ for different angles and spatial resolutions). The former can be seen within 
the flux tube (and thus only at high spatial resolution) and is associated to the lifetime of 
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Line width Mg II 2796.35 Optically thin case Resolution = 0.0 R 
Variation along 0° LOS piane 



Variation aiong 15° LOS piane 



Variation aiong 30° LOS piane 



Variation aiong 45° LOS piane 



Variation aiong 60° LOS piane 



Variation aiong 75° LOS piane 



Variation aiong 90° LOS piane 
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Fig. 7.— Time slice of the line broadening in Mg II k for the case of optically thin intensity 
along a slit perpendicular to the prominence axis, located at the apex. The line broadening 
is calculated from single Gaussian fits and corresponds to the FWHM of the spectral line. 
The spatial resolution is 0.1 R. Each panel corresponds to a different LOS angle. 


KHI vortices, on the order of one period or less. The latter component is associated to the 
heating of the entire flux tube and occurs in a timescale of 10 — 15 min (roughly two periods, 
between times t = 10 and t = 25 min), and a few minutes shorter in the case of the high 
resolution run. 
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Intensity Si IV 1402.77 Optically thin case Resolution = 0.3 R 
Variation along 0° LOS piane 
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Fig. 8.— Time slice of the intensity for Si IV along a slit perpendicular to the prominence 
axis, located at the apex. The spatial resolution of 0.3 R. Each panel corresponds to a 
different LOS angle. 


Although a proper parameter space study is needed for investigating the full effects of 
the initial transverse wave amplitude, results from the previous study ( [Antolin et ah 2014) 
indicate that an increase (decrease) of amplitude (within the linear regime) produces stronger 
(weaker) current sheets and larger (smaller) vortices, due to the increase (decrease) in velocity 
shear. Stronger (weaker) ohmic and viscous dissipation and turbulence is thus obtained. We 
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can then expect that chromospheric plasma within the small-scale vortices composing the 
turbulence will tend to be more or less rapidly heated depending on the amplitude of the 
oscillation. Specifically, smaller amplitudes produces longer lived flux tubes (as observed in 
chromospheric lines). This model can therefore explain the observed tendency of lifetime 
decrease with height, shown in Fig. 2 of Paper 1, since amplitude of the transverse motion 
is observed to increase with height. It is also likely that other longitudinal effects besides 
heating, such as the length of the flux tube or the mass flows, play an important role in the 
morphology and lifetimes of threads (Zhou et al. 2014). 


4.2. Damping of transverse oscillations 


Although signatures of damping can be seen in the threads, constant amplitude oscilla¬ 
tions and even sudden increases can also be seen in movies 3 — 6 of Paper 1. Such features 
may seem at first contradictory with the resonant absorption scenario, since this one indi¬ 
cates that a net transfer of energy should take place from the global kink mode (affecting 
the transverse POS motion) to the local Alfven waves in the resonant layer (affecting the 
azimuthal flow and therefore the LOS velocity component). However, there are a number of 
realistic scenarios that indicate that damping is not a necessary observational characteristic 
of resonant absorption. For instance, a decrease in the density of the flux tube or the back¬ 
ground can lead to an increase of the amplitude of the oscillation, thus counteracting the 
damping from resonant absorption (in other words, an increase in amplitude does not imply 
an increase of energy). Similarly, a decrease of the flux tube’s density or an increase of the 
background’s density increases the damping time. A decrease in density would also lead to a 
decrease of the opacity of the material, making the oscillation hard to follow. Another possi¬ 
bility for non-decreasing kink amplitudes is one in which the oscillation energy is ceaselessly 
injected into a damped harmonic system through a non-resonant continuously operating ex¬ 


ternal force. This has been suggested by Anfinogentov et al. (2013) based on undamped 


transverse oscillations in loops observed with AIA. Such external force could be readily (and 
is most likely) achieved by continuous convective motions at the photosphere. In this sce¬ 
nario, the energy that is continuously dumped into the system is damped by the resonant 
absorption (that we observe) and turned into heat, without actually changing the amplitude 
of the wave significantly. Another possibility for the apparent decay-less oscillations is pro¬ 
vided by the numerical model and suggests that the observed threads may correspond to 
the obtained threads in the simulation, which are the result of KHI vortices. The observed 
POS oscillations would then correspond to those of the KHI vortices, which, as explained in 
section [4^ produce oscillating thread-like structure that damps, as an ensemble, on longer 
timescales. 
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4.3. X-y symmetry 

The characteristic arrow shaped structure seen in the Doppler maps of Fig. implies 
that the resulting X — y profile is generally symmetric, where y denotes the distance across the 
flux tube in the POS. Indeed, by taking a specific time in the figure and plotting the Doppler 
velocity versus perpendicular distance, the resulting profile is symmetric with respect to the 
center of the flux tube. Such symmetry has strong consequences since it implies that in an 
observational setup such as that in Paper 1 (Fig. 5) a coherent signal between the slits is 
expected. To better understand why this is so, we have produced Fig. [^with RH from a 
random snapshot of our simulation (same as in Fig. 7 of Paper 1). Similar to Fig. the 
figure shows that when viewing along an inclined ray that has a significant perpendicular 
component with respect to the direction of oscillation, the Doppler signal is to a large extent 
symmetric across the flux tube (as projected onto the plane of the sky). In the fixed x — y 
frame of reference, in which the kink mode oscillation occurs along the x-direction, the 
y components of the azimuthal flows have different signs (thin white arrows) and become 
zero at the trailing and leading (with respect to the direction of kink oscillation, i.e., the 
x-direction) edge of the flux tube. For our spectral lines (Mg II k and Ca II H), rays 1 
and 3 which cross the edges of the flux tube will go through only optically thin material, 
as depicted in Fig. For ray 3 and basically any viewing angle there is an imbalance of 
red-shifted and blue-shifted material resulting in a blueward Doppler shift. This imbalance 
is caused by the fact that the azimuthal flows are zero at the trailing and leading edges of 
the flux tube. A similar imbalance in ray 1 causes the same type of blueward Doppler shift. 
This explains the general symmetry of the velocity across the flux tube observed in Fig. 5 of 
Paper 1. There are small deviations to this symmetry (e.g., the jagged nature of the X — y 
plot between rays 2 and 3 in the right panel of Fig. |^. These deviations are caused by 
the fine-scale velocity signals associated with the KHI vortices on the boundaries of the flux 
tube. This is similar to what we see in the observations (Fig. 5 in Paper 1). We note that in 
Fig. 11 the optical thickness of the flux tube in Mg II k is automatically taken into account 
and is thus not a free parameter. This is important for the comparison of our model results 
with the observations. 

Only if the tube were completely optically thick, for very small shallow angles and in the 
case of a flux tube whose internal flows were resolved with the spectrometer (case in which 
the flows from the individual KHI vortices can be resolved) the symmetry of Doppler signals 
across the tube could be significantly changed because then only the side of the flux tube that 
faces the observer (and potentially the individual KHI vortices) contributes to the observed 
signal. In addition, the presence of the KHI thickens the layer in which significant azimuthal 
flows occur, thus facilitating a more symmetric Doppler pattern across the flux tube. Both 
of these effects cause significant differences between our model and previous models who 
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considered completely optically thick flux tubes and did not include the KHI (Goossens 


eFahllMil ). These different results clearly show the importance of proper radiative transfer 


modeling in order to forward model numerical results. 



Fig. 9.— Doppler shifts variation across the flux tube in the Mg II k line. \ — y diagrams 
(panel b) produced with RH in the Mg II k line from a slit placed at the center and perpen¬ 
dicular to the tube’s axis for several inclination angles (in negative colors). The diagrams 
are made from the same snapshot as in Fig. 8 of Paper 1 (shown here in panel a, rotated 
for comparison with Fig. [^. The optical thickness values (yellow contours) are reproduced 
from that flgure. The curved white arrows on panel a indicate the direction of the large 
scale azimuthal flow at that instant in time. Three rays denoted as red (1), yellow (2) and 
green (3) denote 3 examples of LOS rays going through the flux tube at an inclination of 
32°. The Doppler shifts at these 3 locations along the slit are marked on panel b in the 
corresponding panel of that angle. Solid and dashed arrows indicate Doppler shift features 
generated by the resonant flow and the KHI vortices, respectively. Except for the individual 
Doppler signals from the KHI vortices (which can have a sporadic different behavior) the 
Doppler signal across the flux tube is generally symmetric with respect to a LOS through 
the center of the flux tube: the leading and trailing edges (rays 1 and 3) exhibit a similar 
Doppler signal. 


4.4. Dynamical coherence in the transverse MHD wave model 


As seen in Fig. and explained in sections |3.4| and |3.5[ the thread-like structure and 
the characteristic arrow (>) shaped structures can be observed for basically any LOS angle 
and optical thickness in the present transverse MHD wave model. Such behavior is obtained 
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Fig. 10.— Time slice of the optically thin intensity (top left panel) and LOS velocity (top 
right panel) in the Mg II k line along a slit perpendicular to the prominence axis, located 
at the apex and for a LOS angle of 45°. The corresponding time slice intensity in the 
Si IV 1402.77 A line is shown in the bottom left panel. We set a spatial resolution of 0.5 R 
for the Mg II k intensity and 1.0 i? for the LOS Mg II k velocity and Si IV intensity, 
simulating, respectively, what SOT and IRIS would observe in the case of an unresolved flux 
tube (here we assume the Ca II H intensity of SOT to be similar to the optically thin Mg II k 
intensity). In the bottom right panel we show the Doppler signals along the green, orange 
and red slit locations, with respective colors (whose colored edges are shown in the top 
panels), centered at a distance of 3 R, 0 R and —3 R from the flux tube’s axis, respectively. 
The black dashed line corresponds to the transverse displacement of the flux tube in the 
POS, calculated from Gaussian fits to the intensity image for each time step. Times without 
Doppler signal correspond to regions for which the signal is too low. 


by the combination of resonant absorption and the presence of KHI, and is predicted in the 
observations by our advanced radiative transfer calculations. As explained in the previous 
section, the obtained arrow shape in the Doppler maps means a \ — y profile that is generally 
symmetric, where y denotes the distance across the flux tube in the POS. This implies that 
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in an observational setup such as that of Fig. 5 in Paper 1 a coherent signal between the 
slits is expected. 


In a model with negligible resonant absorption, and therefore negligible azimuthal flow in 
the resonant layer, we would have the transverse displacement and the dipole-like azimuthal 
flow around the tube. Because of the little transfer of energy from the transverse motion 
into the resonant layer in such case, the plasma with azimuthal flows would mainly be 
external plasma, which we expect to be at a considerably different temperature (probably 
coronal) than the temperature within the flux tube. This would therefore be invisible to 
IRIS. However, as shown by Fig. 5 in Paper 1 the out-of-phase relation between the POS 
displacement and the LOS velocity exists over a significant distance across the threads, 
suggesting that the resonant layer is not bound to the thread but can be detected over 
a larger distance. This result therefore cannot be explained by a transverse MHD model 
without a resonant absorption layer (or one in which the resonant layer is highly confined to 
the boundary layer and is not broadened by any other mechanism). On the other hand, as 
shown by Figs. and (see also Fig. 6 in Paper 1 and Fig. 10), in our numerical model 
combining resonant absorption and KHI the characteristic Doppler signal extends over a 
significant transverse distance with respect to the corresponding chromospheric intensity 
signal from the flux tube, matching the observations. 


Two effects contribute significantly to this spatial broadening. First, the KHI rapidly 
broadens (inwardly) the boundary layer after it sets in, clearly seen in Fig. |^and movie 1. 
The KHI extracts the energy out from the resonant layer and imparts momentum to the 
eddies around the boundary layer. As a result, the plasma layer with azimuthal flow has 
a significant width, which can then be detected with the spectrometer. As explained in 
section |3.5| this spatial broadening is not observed at high resolution in the optically thin 
intensity images in Mg II (with which actually a thinning of the flux tube is observed), but 
in those of Si IV, as shown by Figs. |^and|^ and in the Doppler and line broadening maps, 
as shown by Figs. and 


The second effect is perhaps even more significant and relies on instrumental aspects: 
when calculating intensities and velocities from our numerical model we have taken into 
account the different spatial resolutions of SOT and IRIS for correct comparison with the 
observations in Paper 1 (further explained in section 4.6). The coarser resolution of IRIS 
(0.33" — 0.4") with respect to SOT (0.2") implies that the Doppler signal spreads over a 
wider distance than that delimited by the SOT threads, and corresponds at least to an area 
proportional to that of the IRIS PSF. Since the resolution of IRIS is roughly the same as the 
width of a thread (as shown in Fig. 5 of Paper 1), the signal should be detected at least a full 
diameter away from each side of the thread. This closely matches what we observe, as shown 
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by Fig. 5 in Paper 1, thereby supporting the fact that our Doppler measurements with IRIS 
in Paper 1 actually are set by the resonant absorption region (and not the exterior). 

It is further important to have in mind that Mg II h&k emission exists all around 
the SOT threads, as evidenced by Fig. 1 in Paper 1. This, together with the existence of 
dynamical coherence, may suggest a scenario in which the threads belong to a larger inho¬ 
mogeneous flux tube (which would exist over the distance in which coherence is detected). 
Such a scenario is discussed in section 14.61 


4.5. Phase mixing 


In our simulations two important mechanisms for small-scale generation are present: 
KHI and phase mixing. However, the timescales between KHI and phase mixing are different 
when it comes to produce small scales. Indeed, from one hand, over several periods phase 
mixing generates multiple very small-scale layers of counter-streaming flows located around 


the resonant layer. This also occurs in the case of an initially inhomogeneous loop (Terradas 


et al. 2008c). By itself, phase mixing would then result in a broadening of the line width. 


a decrease of the Doppler signal and a loss of coherence (this is clearly seen after 6 full 


periods in Terradas et al. 2008c). Such behavior is also confirmed by our numerical model. 


Phase mixing produces the out-of-phase behavior rapidly (due to the change of the Alfven 
speed across the boundary layer), while its generation of smaller scales (counter-flows) is 
produced on longer time scales. Accordingly, a decrease of the Doppler signal is clear when 
comparing the end stages with the first stages of the simulation. On the other hand, the 
small scales produced by KHI occur rapidly (within a period) and line broadening from KHI 
is readily produced, as seen in Fig. Importantly, the small-scale vortex flows produced 
by the KHI near the boundary receive momentum from the larger scale azimuthal flow from 
the resonance, whose characteristics therefore dominate the average signal picked up by the 
slits, thus avoiding the loss of Doppler signal. 


4.6. Effect of spatial resolution - KHI vortices as threads? 

The transverse sizes of the observed thread-like structures with SOT vary between 200 
and 400 km approximately. As discussed previously, in phase oscillation of threads, both 
in the POS motion and in the associated LOS velocity, exists (movies 3 — 6 and Fig. 5 in 
Paper 1), suggesting that some dynamic coherence exists in the transverse direction with 
a length scale of a few 1,000 km or so. Assuming a resolving power of 0.1 R for SOT and 
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given that the Ca II H emission is similar to that of Mg II k, but optically thin with little 
opacity (as shown in Fig. 7 of Paper 1), we would expect an SOT time slice to be similar to 
Fig-i As shown by the figure, in this case SOT would be resolving the internal structure 
of the prominence flux tubes. The fact that Mg II h&k emission exists all around the SOT 
threads (as shown by Fig. 1 in Paper 1) further suggests that these structures may belong 
to a larger flux tube. It is therefore reasonable to consider that the observed thread-like 
structure would be the result of KHI vortices. Keeping roughly the ratio between the spatial 
resolutions of SOT and IRIS (a ratio of 2 —3), in this case the spatial resolution of IRIS would 
be 0.3 R. Fig. 6 in Paper 1 would therefore apply to this case, meaning that the observed 
characteristic coherence in the transverse direction and the phase difference between the 
transverse displacement in the POS and the LOS velocity can be explained. 


As shown in Figs. and the KHI vortices are short lived (on the order of a period or 
so). This short lifetime can explain the observed fast disappearance of threads and makes 
the damping hard to observe in the numerical model as well. Notice that the damping of 
the transverse POS motion of the flux tube is only clearly present at the beginning of the 
simulation. This is because of the fast damping of the quasi-mode, which also holds for 
an inhomogeneous multi-threaded loop (Terradas et al. 2008ci Pascoe et al. 2011). The 


presence of KHI further complicates the detection of the overall damping due to the large 
amount of thread-like structures and line-of-sight superposition that occurs (please keep 
in mind that, as suggested by our numerical model the Ca II H line may be significantly 
optically thin and therefore a significant superposition of such threads is expected in the SOT 
images). Furthermore, the KHI is fed by resonant absorption and the KHI vortices are partly 
embedded in the azimuthal (dipole-like) flow around the boundary layer that results from 
the resonance (movie 1). This implies that the oscillations of the KHI threads can increase 
in amplitude and that the damping is much weaker than that of the POS motion of the flux 
tube since the energy is continuously fed in via resonant absorption. This is especially so 
since at the low resolution of the current instruments (compared to the numerical resolution) 
and since they move coherently to a large extent due to the large scale azimuthal flow in 
which they are embedded each thick thread is actually an unresolved ensemble of threads. 
This behavior can be seen in Figs. and for Mg II k and in Fig. |^for the Si IV line, and 
matches well the observations (section |4.2| ) . The damping of the ensemble of threads is only 
clearly observed when comparing the general ensemble of oscillations by the threads at the 
beginning and at the end of the simulation, and depends on the continuous generation of the 
KHI at each oscillation and the cascading of energy towards the lower scales. The former is 
limited by the energy provided by the resonance, and the latter is guaranteed by the high 
Lundquist number in the corona. Note that being able to resolve the KHI threads does not 
mean that we are fully resolving the boundary layer of the flux tube. The KHI threads retain 
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in general the same phase with the POS motion of the flux tube, as shown in Figs. and 
1^ The detected Doppler signal, coming from a coarser instrument than the imaging signal, 
does not correspond to the individual signals of each thread but rather to the average signal 
of the ensemble of threads, as shown in Fig. 6 of Paper 1. 

Despite the many observational facts supporting the interpretation of KHI vortices as 
threads, we must also consider the case in which we are not resolving the oscillating flux tube. 
In this case, we set a spatial resolution of 0.5 R and 1.0 R for SOT and IRIS, respectively. 
As can be seen in the top left panel of Fig. [T^ little internal structure is obtained, which 
suggests that the observed thread-like structure is not related to KHI in this case. However, 
the panels on the right hand side clearly show that the characteristic phase difference is 
maintained. Interestingly, even at this coarse resolution the ensemble of threads can be 
picked up in the hotter Si IV line as a single large thread (shown in the bottom left panel). 
This single ‘thick’ thread produced by the ensemble of KHI threads conserves a very similar 
oscillation phase as the Mg II k intensity oscillation (and is therefore out of phase with 
the LOS Mg II k velocity), shows little damping due to the momentum transfer from the 
azimuthal waves, and thus provides good agreement with the observations. 


5. Conclusions 


In this work we combined high resolution observations with Hinode and IRIS (Paper 1) 
with advanced 3D MHD numerical resolutions coupled with appropriate forward modeling to 
present the observational signatures (imaging and spectroscopic) of transverse MHD waves 
in prominence flux tubes. Resonant absorption rapidly produces characteristic out-of-phase 
behavior between the POS motions (produced by the kink quasi-mode) and the LOS velocity 
(produced by the azimuthal flows generated by the torsional Alfven waves). The azimuthal 
flows from the resonance further produce localized line broadening at the boundaries of the 
flux tube (also through the process of phase mixing). These features would be invisible 
to current imagers and spectrometers if it were not for a secondary phenomenon that can 
accompany transverse MHD waves: KHI instabilities rapidly produced in timescales of one 
period by the velocity shear at the boundaries of the flux tube. A complex interplay between 
resonant absorption and the KHI occurs in which the KHI extracts the energy from the 
resonant layer and dumps it into heat through viscous and ohmic dissipation at the generated 
vortices and current sheets. The vortices, combined with LOS projection effects (especially 
for optically thin emission) lead to thread-like structure within the flux tube with enhanced 
emissivity in chromospheric lines, similar to the coronal case ( Antolin et al.||2014 ). The large 
vortices rapidly degenerate into turbulent flows that propagate inwardly and enlarge the 
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transition boundary of the flux tube between the core and the external medium. A major 
result is that since the KHI vortices receive momentum from the resonance flows the out- 
of-phase behaviour and line broadening from the resonance are transferred to the enlarged 
transition boundary, therefore making these phenomena observationally detectable. We have 
shown that the observational signatures from this model match well the observations. 

In the transverse MHD wave interpretation that we have put forward for explaining the 
observations we have considered two different scenarios depending on whether the flux tube 
in which the wave exists is resolved or not. The first case considers the observed threads as 
independent flux tubes for which we are not able to resolve their internal structure. Each 
thread undergoes a kink oscillation and the POS motion and LOS velocities we observe are 
associated to each kink separately. A sketch of this case is shown in Fig. 8 of Paper 1, 
panel A. In this work we have considered this case modeling a single flux tube and shown 
that the observed characteristic coherence over a significant transverse scale and the out-of- 
phase relation between the transverse displacement in the POS and the LOS velocity can be 
explained. The obtained timescales for the heating of the flux tube in our model match well 
the variability in the emission that can be associated to heating in the observations. We have 
shown that these features are numerically robust. Indeed, when increasing the resolution 
and Lundquist numbers an increase of small-scale vortices and current sheets is obtained, 
which enhances the heating due to KHI but keeps the large-scale dynamics unaltered. The 
results are also observationally robust, being to a large extent independent of the optical 
thickness of the material and the LOS angle (except in the completely optically thick case 
and very shallow angles). The observed dynamic coherence and out-of-phase relation are 
an imprint of resonant absorption and depend on the KHI only for widespread heating and 
thus detectability. At high spatial resolution and in the case of a well-defined flux tube 
the turbulence generated by the KHI is reflected in an increase of line broadening towards 
the boundary, especially when observing along the direction of oscillation. The detection 
of those dynamic features over a wide range of LOS angles and spatial resolutions is thus 
greatly facilitated by the presence of this instability. 

The second case corresponds to a scenario in which the internal structure of the flux 
tube is resolved. In the previous section we have suggested that in this case the observed 
threads could be the KHI vortices, and shown that the expected signatures match well 
with the observed SOT threads. A confirmation or rejection of this scenario may only be 
achieved with higher instrumental spatial resolution. Another suggested possibility is that 
of an inhomogeneous flux tube in which regions of high density within the flux tube host 
the observed threads. This case is pictured in Fig. 8 of Paper 1, panel B, and we can refer 
to it as the multiple kink wave case. Each high density region can have its own resonant 
layer, leading to azimuthal flows around each thread. Numerical simulations of this case 
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but in a fully coronal setup with a multi-stranded loop with an irregular cross-section clearly 
indicate that the ensemble of strands (threads in our case) evolve coherently, with a dominant 
frequency everywhere within the structure in a similar way as for the initially homogeneous 


loop case (Terradas et ah 2008c; Pascoe et ah 2011). In other words, a large azimuthal flow 


also develops around the combined strand (thread) system. These results strongly support 
the presence of a larger scale azimuthal flow around the threads with the same kind of out-of- 
phase behavior between the POS motion and the LOS velocity. Our results for the initially 
homogeneous flux tube can therefore also be applied to this initially inhomogeneous flux 
tube case. 

Lastly, we have shown that the heating produced by the KHI affects a large area around 
the (smaller) prominence core, bringing the temperature of the plasma in this layer from 


chromospheric to transition region values, an effect also reported in the coronal case (Antolin 


et al.||2014 ). This signature is in turn picked up at any LOS angle by TR lines such as Si IV, 
leading to a gradual fading of the threads out of the chromospheric pass bands (except for 
the prominence core) and appearance in the TR pass bands. This effect is accompanied by 
an apparent thinning of the flux tube, especially in optically thin emission in chromospheric 
lines such as Ca II H. 

All the mechanisms presented here for prominence plasmas extrapolate also to the coro¬ 


nal scenario (of which only the imaging aspects are shown in Antolin et ah 2014), and may 
also extrapolate to spicules, in which multi-threaded structure, transverse swaying and tor¬ 


sional motions have recently been reported (De Pontieu et ah 2014a; Rouppe van der Voort 


et ah 2015; Skogsrud et ah 2014). The combination of resonant absorption and KHI may 


therefore play an important role in the chromospheric and coronal morphology and heating. 
The assessment of this mechanism as such is the subject of future work. 
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